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Abstract. The W Mass and QCD Working Group discussed a wide variety of 
topics relating to present and future measurements of Mw at LEP 2, including 
QCD backgrounds to W~^W~ production. Particular attention was focused 
on experimental issues concerning the direct reconstruction and threshold mass 
measurements, and on theoretical and experimental issues concerning the four jet 
final state. This report summarises the main conclusions. 



1. Introductionf 

The 'W Mass and QCD Working Group' addressed a variety of topical questions during 
the Workshop. The format varied from formal seminar presentations to informal 
discussions, with a total of 28 people contributing. In this article we summarise the 
outcome of the discussions, including in particular new results obtained both during 
and after the meeting. We do not attempt to review the status of the various physics 
topics prior to the meeting, as this was very well covered in the plenary talks by J C 
Thompson B R Webber § and G Cowan [|. 

For most of the time the Working Group separated into two partially overlapping 
subgroups. The first focused on theoretical and experimental issues concerning various 

f Unless otherwise stated, the sections have been prepared by the convenors. 
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aspects of the final state in W^W~ and qcd four jet production, in particular colour 
reconnection, Bose-Einstein correlations, and the accuracy of current qcd models for 
the four jet final state. The second subgroup was concerned with mainly experimental 
issues concerning the direct reconstruction and threshold cross section methods for 
measuring Mw at lep 2. In addition, the subgroup updated the expected precision of 
the two methods based on experience with the two methods to date. 

The report is organised as follows. The work of the two subgroups is described in 
Sections ^ and ^. Each section contains a general overview followed by individual 
contributions as subsections. The overall conclusions of the Working Group are 
presented in Section ^. 



2. Aspects of the hadronic final state in W^W production 

2.1. Experimental aspects of colour reconnection § 

Colour reconnection (also referred to as 'rearrangement' or 'recoupling') in W+W" 
decays has been the subject of many studies (e.g. @, ||, ^) and at present there is 
agreement that observable effects of interference between the colour singlets in the 
perturbative phase are expected to be small. In contrast, significant interference 
in the hadronisation process appears a viable prospect but, with our current lack 
of knowledge of non-perturbative qcd, such interference can only be estimated 
in the context of specific models [f[ @, |[ |lO[ In the studies described 

below, experimentally accessible features of these models|| are investigated, paying 
particular attention to the bias introduced to a typical measurement of Afyv by direct 
reconstruction of the decay products. 

Throughout this section reconnection effects were studied using: pythia 5.7 
type I and type II superconductor models (with the string overlap integral in the type 
I case characterised by p = 0.9) ^ ||]; ariadne 4.08 allowing reconnection between 
the two W bosons; and herwig 5.9, in both its default reconnection model and also 
a 'colour octet' variant in which merging of partons to form clusters was performed 
on a nearest neighbour basis^. In all cases, the tuning of the models was as used in 
reference [ p2| . 

2.1.1. Inclusive charged multiplicity It has been suggested Q that simple 
observable quantities such as the charged multiplicity in restricted rapidity intervals 
may be sensitive to the effects of colour reconnection. More recently |11 it was 
suggested that the effect on the inclusive charged multiplicity itself may be larger 
than previously considered and that the mean hadronic multiplicity in W"'"W~^qqqq 
events, (N^i^), may be as much as 10% smaller than twice the hadronic multiplicity in 
W+W^— >qq£Ff events, {N^^^'^)- It was also reported during this workshop that the 



effects of Bose-Einstein correlations may increase 

(^ch> by - 3-10% (see section 
The shifts in {N^l^) at the hadron level predicted by the models studied thus far 
are given in table 0, where A{N^^) is defined as the change in mean multiplicity 
relative to the 'no reconnection' scenario of each model. From these, it is clear that 
the multiplicities themselves and also the magnitude and sign of the predicted shifts 
are model dependent. 

§ Prepared by M F Watson, N K Watson 

II In studying these models, no retuning was performed when reconnection was enabled. 
^ This was suggested by B R Webber, as a partial emulation of the model of reference 
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Table 1. Mean charged multiplicities, (1^^^), and predicted shifts for various models 



model 








MN^S) {%) 


PYTHIA 


normal 




38.64 






type I 




38.21 


-l.litO.l 




type II 




38.39 


-0.7±0.1 


HERWIG 


normal 




37.07 






reconnected (P = 




37.25 


+0.5±0.1 




reconnected (P = 


1) 


38.38 


+3.5±0.1 


ARIADNE 


normal 




38.14 






reconnected 




37.07 


-2.8±0.1 



In this study, the precision with which such tests may be performed is quantified. 
As a starting point for such tests, it was first verified that in the absence of 
reconnection effects (N^^) = 2{N^^^'^) in the models pythia and herwig. The 
statistical uncertainty of this test was 0(0.1%). Next, samples of 10^ herwig and 
PYTHIA W+W^ events were generated at y/s — 171 GeV with a full simulation of the 
OPAL detector, and realistic event selections were applied for both W+W^^qqqq and 
W+W~— >qq£Ff (i = e, n and t). The efficiency in each case was ~80%, while the 
purity is ~ 80% for W+W'^qqqq and - 88% for the W+W~->qq£l7£ channel. 

The resulting (uncorrected) charged multiplicity distributions for the hadronic and 
semi-leptonic channels are shown in Figs. 0(a) and^b), respectively. The simulated 
data correspond to an integrated luminosity of 10 pb~^ at ^/s = 171 GeV, i.e. that 
delivered by lep in 1996. In both distributions, the expected background is shown as 
a hatched histogram. The significant level of Z''/7 ^ qq background is apparent in 
the fully hadronic channel. 

To extract the mean charged multiplicity at the hadron level at a fixed centre- 
of-mass energy from such distributions, one can apply a simple correction, based on 
Monte Carlo, to the observed mean value, after subtracting the expected background 
contribution. An alternative is to carry out a matrix-based unfolding procedure using 
the cvcnt-by-event correlation between the charged multiplicity at the hadron level 
and that observed in the detector after all analysis cuts have been performed. A 
separate correction for the effects of initial state radiation are necessary in this latter 
case. A third alternative is to integrate the fragmentation function but this is not 
discussed here. 

Based on the the simulated data in Fig. |l|(a) and (b), the expected statistical 
uncertainty on the difference (N*^) — 2{N^^^'^) for an integrated luminosity of 10 pb"""^ 
is 2.2 units, or 5.7% on (N^^). The evolution of the precision of such difference 
measurements with more data is estimated using the following assumptions. Firstly, 
the distributions of N^^ and N^^^'^ are seen to be relatively insensitive to changes in 
centre-of-mass energy once away from the threshold region, as illustrated by the energy 
dependence of (N^^) in Fig. |l](c). Therefore both the mean and the corresponding 
rms are assumed constant at their 184 GeV values. Secondly, above ^/s = 184 GeV 
the W+W" production cross-section is predicted to vary by less than 10% in the 
region up to < 200 GeV, and so a constant cross-section of 16 pb is assumed. 
Thirdly, it is assumed that the selection efficiency at 171 GeV may be maintained 
at higher energies. The expected background cross-section is not important as it is 
subtracted in performing the measurement. Given these assumptions, the dependence 
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Figure 1. Inclusive charged multiplicity distributions with 10 pb ^ of fully 
simulated data, with background indicated hatched, at ^ = 171 GeV for (a) 
W+W~^qqqq, and (b) W+W"^qq£i7f events, (c) Variation of (N*^) with y/s. 

(d) Luminosity dependence of the statistical uncertainty of {N^^) — 2{N^^^'^) (units 
of multiplicity). 



of the expected statistical error on the difference, 6{{N^^) - 2{N^^'"')), is shown as a 
function of integrated luminosity in Fig. |^(d) . 

Typically in such multiplicity determinations, systematic effects become significant 
below a statistical precision of 0.5 units of multiplicity. Uncertainty in the modelling 
of 4-jet like Z°/7 — > qq background with parton shower Monte Carlos in the fully 
hadronic channel may become a significant systematic. 



2.1.2. Event shapes Global event shape variables have been considered in earlier 
studies as potential signatures for reconnection |||, |ll) . In most studies the predicted 
effects on such observables induced by reconnection has been sufficiently small that 
detection would be marginal, even with an integrated luminosity of 500 pb^^. 

The choice of a 'no reconnection' reference sample with which to compare data 
deserves some thought. In trying to find sensitive observables, using the models 
alone is ideal. However, once possible signatures have been developed, and one 
starts to search for effects in data, it will be invaluable to have a well defined 'no 
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Figure 2. (a) Effect of typical experimental selection on thrust distribution, and 
(b) hadron level rapidity distribution in ARIADNE for p < 1 GeV. 



reconnection' reference sample in data to reduce model and tuning dependence, lep i 
data provide a high statistics reference, but additional assumptions are necessary in 
either extrapolating energy scales, or in combining pairs of — > qq to emulate 
W^W^— >qqqq events without reconnection. It is also necessary to assume that data 
recorded and processed by the detectors before 1996 can be directly compared with 
those recorded near the end of the lep 2 programme. For some signatures, the ideal 
reference data are W+W~— >qq£F£ events. However, this sample has only limited size 
and the comparison may require the association of pairs of jets with Ws in the fully 
hadronic channel, a procedure which experimentally introduces more uncertainty. In 
the following, all changes are relative to the 'no reconnection' version of each Monte 
Carlo model and all samples are W+W~^qqqq. 

This study compares the differences in the rapidity distribution of charged 
particles, dNch/dy^ relative to the thrust axis of each event, in the central region, 
|y| < 0.5 and for all y, as suggested in j|, ||, Q]. As the effects are expected to be 
more pronounced for softer particles, the distribution is studied for three momentum 
ranges, p < 0.5 GeV, p < 1 GeV and all momenta. It has been suggested 0, |ll| that 
reconnection effects may be more pronounced in specific topologies where the quarks 
from different Ws are close to one another, therefore events are also studied for all 
thrust values and for T > 0.76. One aspect not considered in previous studies has 
been the effect of applying a realistic event selection, which is necessary to reduce 
the large background {<7{7P qq) ~ 20cr(W"'"W~— +qqqq)). As this is dominated 
by two-jet like events, the efficiency for selecting W"*"W^— ^qqqq events in a similar 
configuration is relatively small, as illustrated in Fig. ^(a); ^ 38% of W+W^^qqqq 
events selected satisfy T > 0.76, falling to < 0.05% for T > 0.92. 

In ||ll[], the rapidity was studied relative to the axis bisecting the two di-jet axes, 
as a function of the angle separating these axes. Experimentally, without any reliable 
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charge identification algorithm to separate quarks from anti-quarks, the specific angle 
proposed in must at best be folded in experimental analyses, and also requires 
pairing of jets into Ws. While the rehability of associating the 'correct' jets together is 
possible with moderate efficiency using kinematic fits, selecting high thrust events was 
used in the current studies for expediency and simplicity. As the shifts in Mw expected 
are modest compared to the experimental mass resolution on an event-by-event basis, 
it is worth considering the use of kinematic fits in which our current knowledge of 
Mw is applied as a constraint, in a similar way to that used by experimental TGC 
analyses. 

Hadronic events were generated using the models pythia, herwig and ariadne, 
with and without a simulation of the opal detector, and dNch/dy studied within the 
ranges of y, p and T described above. A smearing simulation of the opal detector, 
which is reliable for studies in the W+W^— s-qqqq channel and necessary to achieve 
the relatively high statistics required, was used herein and also to estimate shifts in 
Mw- 

As an example of how the differences may be concentrated in restricted rapidity 
intervals. Fig. ||(b) shows the dNch/dy distribution for p < 1 GeV in ariadne, for 
events with and without reconnection. Changes in charged multiplicity, A(A^^j^), 
within given p and y intervals are summarised in Fig. ^(a) for each of the models 
introduced in table ||, without detector simulation. The left (right) hand side of the 
figure shows the percentage change in (N^^) for the three momentum ranges considered 
for all y {\y\ < 0.5). The leftmost points in this figure correspond to the results of 
table |l]. Fig. |(b) gives the analogous results for T > 0.76. For illustration, statistical 
errors corresponding to an integrated luminosity of 500 pb~^ are given for the 'herwig 
colour octet' model. 

It is seen that in all models the magnitude of the change increases when only low 
momentum particles are considered. Applying a thrust cut such as T > 0.76 rejects 
^ 40% of events and may change (N*^) by up to two units, but differences relative 
to the 'no reconnection' scenarios are essentially unchanged, therefore the sensitivity 
is reduced. The predicted maximum statistical significance of A{N^^), as well as its 
sign, depends strongly on the model, varying from ^ 6a for ariadne and the herwig 
'colour octet' model, ~ 3.5cr for pythia type I, ^ 2ct for pythia type II, down to 
^ O.SfT for HERWIG. The point of maximal sensitivity is indicated (square markers) 
for each model in the figure. Similar trends were observed in studies with detector 
simulation but typically A{N^^) was found to be ~ 50% smaller. 

It may be possible to increase the sensitivity to reconnection effects using charged 
multiplicity based methods, by considering particle distributions relative to the 
W+W~ decay axis, as reconstructed using kinematic fits. In an alternative 
multiplicity signature ('interjet multiplicity') was introduced, having similar sensitivity 
to integrating dNch/dy for |y| < 0.5. This interjet multiplicity was similar in idea 
to methods normally used to quantify the 'string effect' in 3-jet e+e" events. It 
was suggested that this be studied further, using the shape of the particle density 
distribution as a function of the angular separation between jet pairs, rather than 
restricting the study to the integrated particle density in the fixed angular regions. 
However, the 4-jet case is somewhat more complex than the familiar 3-jet case, being 
non-planar, and so this was not pursued during the workshop. 

2.1.3. Shifts in Extracting Mw from the decay products observed by 

experiments is non-trivial, requiring much attention to bias induced from effects such 
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Figure 3. Fractional change in charged multiplicity as a function of maximum 
particle momentum, in two rapidity regions, for (a) all T, and (b) T > 0.76. See text 
for details. 



as initial state radiation, detector calibration, imperfect modelling of the underlying 

physics processes and of the apparatus. In comparison to this, estimating a shift 
which could result from the effects of reconnection phenomena is straightforward, as 
the value of interest is the relative shift between Mw determined in two different 
scenarios of the same model. The absolute value of "Mw" obtained is not central to 
these studies. However, there are still many uncertainties inherent in such studies, 
such as sensitivity of the method used to extract Mw to changes in ^/s, to tuning of 
the Monte Carlo models (e.g. virtuality cut-offs in the parton shower development), 
to treatment of combinatorial background and ambiguous jet-jet combinations, and 
the range over which fitting is performed to name but a few. 

In these studies, the method used to extract Mw followed closely that used 
by OPAL for its preliminary Mw results using 172 GeV data. In this, events with 
detector simulation are first selected using the same procedure as noted earlier. Four 
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Figure 4. Typical mass distribution with fit results with detector simulation and 
event selection. 



jets are formed using the k± jet finder, corrected for double counting of energy 
within the apparatus, and a parametrisation of the errors on the measured jet 4- 
momenta is carried out. A five-constraint kinematic fit, in which the jet-jet pair 
masses are constrained to be equal, is performed for each of the three possible jet-jet 
pairings, event by event. A mass distribution is constructed using the mass from the 
combination having the highest probability from the kinematic fit in each event if this 
has probability greater than 1%. A second entry is also admitted if the second most 
probable fit result has probability greater than 1% and within a factor of three of 
the highest probability combination. The aim of this is to include additional mass 
information for events in which the most probable fit combination is incorrect. In 
such events, these two masses are essentially uncorrclated. A typical mass distribution 
formed by this procedure is given in Fig. ^. 

This method was applied to simulated events from each of the models in turn, 
and the shifts obtained are summarised in table |^, where uncertainties on these shifts 
are statistical. The ariadne model predicts a modest shift in mass of approximately 
50 MeV. No significant shift is predicted by the models pythia and herwig. In an 
earlier study, performed in a similar way, significant shifts were determined The 
PYTHIA and ARIADNE models considered in the present study were also included in Q , 
albeit with different model dependent parameters and looser event selection criteria. 

One quite plausible explanation proposed was that the difference was due to the 
significantly more stringent event selection currently used. It has been shown that the 
current selection preferentially rejects events having two-jet like characteristics, which 
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Table 2. Table of shifts in Afw for each model. 

{AMw> (MeV) 



model 






selected events [e ~ 80%) 


all events 


PYTHIA 


type I 




+18±11 


+11±11 




type II 




-13±11 


-19±11 


HERWIG 


reconnected (P = 


i) 


-16±16 


-19±16 




reconnected (P = 


1) 


+13±15 


+8±14 


ARIADNE 


reconnected 




+51±16 


+51±15 



is where reconnection effects may be expected to be prevalent. The rejection of these 
events does not appear to be the reason for small mass shifts, as a similarly small 
effect is observed when all events are selected, as seen in table |^. 

Many possible sources for the difference were investigated in the context of the 
PYTHIA models. Neither changes in the tuning of pythia/jetset by OPALf to improve 
the description of lep i data, nor the different centre-of-mass energy (i/s = 175 GeV 
in Q) were found to be significant. The current analysis procedure is slightly different 
to that in However, significant shifts are still found when the current procedure 
is applied to the same simulated events used in Conversely, applying the former 
procedure of |^ to the samples herein does not induce a significantly larger mass shift. 

One apparently significant effect was found to be the choice of mass assigned to 
jets in performing kinematic fits. As discussed in Q, this choice is not unique. In 
the analysis of ||^, the hadronic jets were assumed massless whereas in the current 
studies, the measured jet mass was used. Re-analysing the same simulated events 
of ^ but assigning measured masses to the jets reduces the mass shifts estimated, 
e.g. shifts quoted in |] of 130 ± 40 MeV (type I) and 50 ± 40 MeV (type II) become 
70 ± 40 MeV and 30 ± 40 MeV, respectively. For comparison, a sample of 200 000 
fully hadronic type I events were generated at ^/s = 175 GeV using identical model 
parameters and program versions, and analysed using the procedure of also using 
measured jet masses. This yielded an estimated shift of 46 ± 16 MeV. It should be 
noted that fluctuations due to finite Monte Carlo statistics have to be considered when 
comparing with the results of |^ , in which samples sizes for the analogous studies were 
50 000 events. 

Comparing the results for mass shifts in table ^ with multiplicity shifts in table |l|, 
it can be seen that any relationship between them is model dependent. Furthermore, 
relatively large shifts in the charged multiplicity do not necessarily lead to a significant 
shift in Mw- 

2.1.4- Future The future for experimental studies of colour reconnection is quite 
open. There is clear model dependence in signatures and mass shifts may be smaller 
than earlier proposed ||^, although there are other models available ||l^, |l^ which 
were not tested in this study from which different conclusions may be drawn. A 
necessary condition for a model to be taken seriously is that it should describe the 
data, therefore tuning of models has to be addressed. With the current statistical 
precision of lep 2 data, none of the models has been put to a stringent test. The effect 
of background cannot be ignored in the W''"W~^qqqq channel as it proves difficult 



t Among these, the cut-off parameter, Qg, was increased from 1.0 GeV in the similar investigation 
of |||, to 1.9 GeV. 
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to remove. More sophisticated selections may be developed, but typically these make 
use of non-trivial correlations between observables, which may be poorly described by 
the models. A particular concern is the description of parton shower Monte Carlos 
to describe the hard, 4-jet like background which is selected. The remaining point of 
note is that given the model dependence inherent to such studies, it is most important 
to develop signatures which can be tested taking the 'no reconnection' scenario from 
data themselves. 

2.2. Bose-Einstein correlations^ 

As discussed in [|j , studies of the influence of Bose-Einstein correlations on the W mass 
measurement were carried out for the cern lep 2 workshop H, Q , which suggested 
that there could be sizeable shifts in the W mass, AAfw ©(+100 MeV). These 
studies were based on the LUBOEI algorithm, which shifts particle momenta after 
generation of the event, and thus requires rescaling. In this report of the present 
workshop, we outline some more recent results using event weighting schemes. 

2.2.1. Event weighting schemes for Bose-Einstein effects The Bose-Einstein effect 
corresponds to an enhancement in the production probability for identical bosons 
to be emitted with small relative momenta, as compared to non-identical particles 
under otherwise similar conditions. Assuming a spherical space-time distribution of 
the particle source, the correlation function takes the form: 



where Q is the four-momentum difference, = —{pi —^2)^1 a-nd p is the absolute 
square of the Fourier transform of the particle emitting source density, with the 
normalization condition p{0) — 1. The incoherence parameter A takes into account 
the fact that, for various reasons, the strength of the correlations can be reduced. 
Often a Gaussian model is assumed for the source density, which leads to 



where R is the source radius. 

In [|l3j the problem of Bose-Einstein correlations was addressed using an approach 
based on assigning weights to the simulated events according to the momentum 
distributions of final state bosons. In this global event weight scheme a shift 
in the reconstructed W mass distribution may arise if the average event weight 
depends on Mw- The use of global event weights is complementary to the local 
reweighting scheme of ref. |Q in the sense that as opposed to the latter, in the former 
the kinematical properties of the events are preserved, while all probabilities and 
multiplicities may change. The method arises very naturally in a quantum mechanical 
approach, where the weight can be constructed as the ratio of the square of the 
symmetrized multiparticle amplitude to the square of the non-symmetrized amplitude 
corresponding to the emission of distinguishable particles. However, the use of global 
event weights leads to a number of conceptual and computational difficulties, which 
must be overcome before any quantitative conclusions can be drawn. 



C{Q) = 1 + Xp{Q) 



(1) 



p{Q) = exp(-i?2Q2) 



(2) 



§ Prepared by V Kartvelishvili, D R Ward 
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One way of forming the weight is to take a product of enhancements C'{Q) for all 
pairs of identical bosons in the event: 

Vi = n ^^(Qn.J (3) 

il ,12 

For high multiplicity events this weight can become extremely large, so that a few such 
events dominate the weighted distributions and lead to unrealistic results. The event 
weights therefore have to be regularized in some way. In order to keep the statistical 
error at a reasonably low level, events with very high weights (higher than some Knax) 
were discarded. The resulting dependence upon Knax was analysed and the results 
were extrapolated to T^nax oo. 

One can also rescale the weight of the event using a single constant Wq : 

V2 = Vi/W^ (4) 

where Wq is slightly larger than 1, and n is the number of pairs in the event (i.e. 
the number of terms in the product in (^). However, for Vmax fixed and reasonable 
values of Wq the scheme gave rise to numerical difficulties, stemming from the fact 
that increasing brings in more and more events from the high weight tail of Vi, 
which leads to large fiuctuations. The results for the shifts in multiplicity and Mw 
using V2 were roughly consistent with those found for Vi . 

A general problem with any weighting method is, that like the local reweighting 
scheme of ref [|lj], they introduce artificial correlations also between non-identical 
particles. In order to counteract this, one can rescale Vi using the weight calculated 
with non- identical pairs: 

V, = V,/V^'^ (5) 

where Vo is the weight calculated according to (||) but for non-identical bosons in 
the same event, while n and m are the numbers of identical and non-identical pairs, 
respectively. This also leads to a better numerical behaviour, as illustrated in Fig. ^, 
which shows the distributions of Vi and V3 for simulated W"'"W~ events at ^/s — 175 
GeV. The high weight tail is much less pronounced for V3 than for Vi . Both fall off to 
a good approximation as inverse powers, with exponents —2.6 and —1.4, respectively, 
which makes it plausible that the sum of all weights converges. 

A different method of constructing the event weight, which is closer to a 
full quantum mechanical treatment, starts from the introduction of a symmetric 
amplitude, which has n\ terms This leads to a weight: 

A'=/V(Qin)p(Q2.j---p(g™j, (6) 

permutations 

where k is the number of times when the first and second indices differ. For (5 = 0, 
A = 1 and n identical particles, equ. (|^) gives a weig ht of 2"("-i)/2, while equ. 
(^ results in the correct value n\. However, for typical hadronic configurations this 
difference is much smaller, and (^) is rarely used because of computational difficulties: 
by limiting the number of identical particles of each type to 8, one loses about 30% of 
events at the Z peak and about 50% of events in W pair production. 

In a recent attempt to overcome these difficulties, in p^ , the n! permutations 
in (^) were divided into sub-classes where exactly k particles have been exchanged, 
and for low energy hadronic collisions (y^s — 10 — 30 GeV) BE effects were shown to 
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10 20 30 40 50 60 70 80 
Weight 

Figure 5. Distributions of the weights Vi (dashed Hne), V3 (full line) and V5 (dotted 
line) for W+W" events at 175 GeV. 



saturate already at fc = 5. However, possible reasons for this improvement are that 
the multiplicities in hadronic collisions in the above energy range are much lower than 
at LEP, and that the source radius used in ||16| w as i? = 1 fm as opposed to 0.5 fm 
measured at the Z peak at lep i and used in |13|. The larger radius corresponds to 
a narrower distribution in Q-space, so that the peak in the Q-distribution of pions 
does not contribute any longer and the number of significantly contributing pairs is 
dramatically reduced. A subroutine LWBOEI for calculating event weights using this 
approach can be incorporated into jetset/pythia, but it is not clear yet how well 
this approach works for e+e~ collisions at lep energies. 

Similarly, another simplified version of (^) was studied in |l^ to assess BE effects 
on W mass measurements at lep. Here, particles were divided into "clusters" of 
neighbours in Q-space, and simple formulae were derived from under certain 
assumptions (see for details). Detector effects and reconstruction procedures 
have also been included. The authors of also use the source radius of 1 fm, and 
do not see any W mass shift due to BE correlations at the level of their statistical 
precision, concluding that BE effect has a negligible influence, below 30 MeV, on the 
reconstruction of W mass. 

All event weights considered above were based on the Gaussian parametrization 
of the particle emitting source. This implies that Vi > 1 for all values of Q. In 



13 



addition to the above weights, a different pair weight was also studied in [|T3|, inspired 
by recent theoretical studies ||l^. Here p in is not required to be always positive: 

cos(eQfl) 

For ^ close to 1, this is very close to apart from becoming slightly negative at 
large Q. The corresponding weight V5 was built in analogy with (^), but with the 
Gaussian (^) replaced by (|^) (the dotted hue in Fig. ||). ^ = 1.15 was found to lead 
to a good overall description. Due to the better numerical behaviour of this weight 
function (the exponent in the power fit is —2.4), it was possible to apply V5 without 
further rescaling. 

2.2.2. Influence of event weighting on Z properties In order to check the self- 
consistency and inherent systematic errors of the method the same weighting 
procedure was applied to the well studied process of Z hadronic decays. For weight 
calculation, A = 1 was used for pions and kaons originating from sources with decay 
lengths CT < 10 fm, and A = otherwise. The source radius R was taken to be equal 
to 0.5 fm everywhere. In Z decays at lep 1 one observes A « 0.3 if all particles are 
considered, 0.4 if only pions are taken into account and 1.0 for directly produced pions, 



while R w 0.5 fm EG, 21 



Various measurable properties of the Z will be affected to different extents, if 
one introduces event weights into the simulation of its hadronic decays. Since the 
partonic states before hadronization are known to be well described by perturbative 
calculations, which do not take into account Bose-Einstein correlations, uncritical 
application of event weights may lead to large inconsistencies with e.g. measured 
branching ratios and relative frequencies of jet multiplicities etc. In order to see how 
serious these effects are and to judge what consequences this has for the analysis of the 
W+W^ events, the precise experimental data from Z decays can be used to check the 
event weighting schemes of Bose-Einstein correlations. Samples of 100000 hadronic 
events at ^/s = Mz and Mz ± 2 GeV were simulated and the weighting schemes 
described above were applied. Table || presents the changes in the charged particle 
multiplicity, in the apparent Z peak position in hadronic vs leptonic decay modes, in 
the branching fractions for charm and beauty decays (i?c and Rb) and in the ratio of 
three- to two-jet events with and without event weighting. 



Table 3. Differences in charged multiplicity, apparent peak mass of Z, branching 
fractions and three-to-two jet event ratio, between weighted and non-weighted events, 
for various weighting systems described in the text. 



Vi V3 Vs 

A(nch> 3.7 ±0.5 1.3 ±0.2 1.8 ±0.2 

AMzo,MeV 8±3 0±3 1±4 

ARc,% -3 ±2 -2 ±2 0±2 

Aflt,, % -26 ±3 -11 ±2 -5 ±2 

A 3jet/2jet, % 80 ± 20 20 ± 5 20 ± 5 



This analysis resulted in the following: 

• The average charged multiplicity has changed. The weight Vi, which was not 
rescaled, leads to the largest increase when compared to the unweighted results, 
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while both and V5 give a smaller increase around 1.5. In all these cases, 
the change can be accommodated by retuning the parameters in the simulating 
program. 

• In principle, event weighting can result in a shift of the apparent Z mass peak. 
However, only Vi yielded a shift of a few MeV, while for V3 and V5 the shift 
is essentially zero. We have not found any significant change of the apparent Z 
width. 

• The pattern of heavy and light quark fragmentation is rather different. Heavy 
quarks produce significantly fewer pairs with small Q, and all BE effects in this 
approach are less pronounced for heavy quarks. Heavy quark events thus obtain 
smaller average weights, which result in changes shown in Table |[ Note that the 
effect for c-quarks is diluted because 6-quark events reduce the overall average 
weight. In order to exclude this artificial flavour dependence in W decays, the 
weighting and rescaling was performed separately for the different decay modes 
of the Ws. 

• The weighting resulted in a substantial increase of jet activity, as measured by 
the three to two jet event ratio. This is however difficult to quantify because 
of its dependence upon the jet finding algorithm and its parameters. The 
numbers shown in Table || were obtained using LUCLUS with default parameters 
(enjoin = 2.5 GeV), corresponding to fairly narrow jets. The effect decreases for 
broader jets and in any case is much less pronounced in W~''W~ production, so 
no attempt was made to correct for it. 

The reproduced correlation functions for the three weighting schemes, Vi, V3 and 
V5 are shown in Fig. ||. Also shown are fits to the form 

iV(l+/3Q)(l + Aexp(-g2i?2)) (8) 

which is often used to parametrize the experimentally observed correlation function 
in Z decays |20, 2^. The dashed line in each figure represents the result of a fit to the 



correlation function of all particles observed in real data from hadronic Z decays ||21[ 
Vi does not reproduce the correlation function well, giving too small values for 
both the incoherence parameter A and the radius R, while V3 and V5 both give very 
reasonable descriptions. 

Hence one concludes that, provided that the different quark final states (and 
possibly the final states with different number of jets) are treated separately, 
application of the global event weighting technique with rescaling of the weight V3, or 
V5 built using the pair weight is not inconsistent with lep 1 data at the Z, whereas 
the direct application of the product of Gaussian pair weights Vi should be treated 
with more care. 



2.2.3. W pair production pythia 5.7 |2^ was used to simulate the process e+e~ 
W'^W~-^-qqqq, and the weighting schemes described above were applied to simulate 
BE correlations. A basic assumption here is that hadronic W and Z boson decays are 
sufficiently similar, so that by using the tuning of the Monte-Carlo model parameters 
that reproduces the experimental data from Z decays at lep, Bose-Einstein effects 
in single W decays are already effectively taken into account in properties such 
as multiplicities and single particle momentum spectra. Only correlations between 
identical bosons originating from different Ws were included, since BE correlations 
within a single W cannot lead to any change in the W mass compared to the 
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Figure 6. Reproduced correlation functions in Z events using the weights Vi, V3 
and V5 . The dashed hnes show the result of a fit to real data from hadronic Z decays 

El. 



semileptonic channel e+e" W"'"W~^qq£F^. Measurement effects like experimental 
resolution and acceptance, reconstruction method etc. were not taken into account. 

The information from the Monte-Carlo was used to assign each final particle to 
the W"*" or the W~ , as in Q . Ws with mass values in the interval 70 GeV < Mw < 
90 GeV were studied at 175 and 192 GeV to assess the energy dependence. At each 
energy, 10^ events were generated, which is about an order of magnitude higher than 
the expected statistics of all four lep experiments combined at 500 pb~^ integrated 
luminosity per experiment. In general, one expects that BE-induced effects in W+W~ 
production should die out at high energies, as the overlap between the two W decay 
volumes decreases. This requires much higher energies than will become available at 
LEP 2, however, and it is likely that the effect will increase with energy in the lep 2 
range |l|. 

The mass distribution of W bosons was built with and without event weighting 
for each of the weights used, and the differences were calculated in the average 
charged multiplicity rich, the mean W mass, Af|^'^'*", averaged over the whole interval 
70 GeV < Mw < 90 GeV, and a fitted M^. The fit was performed using a relativistic 
Breit-Wigner shape with an s-dependent width, in the interval 80.25 ± S GeV, with 
5 — 2 GeV. The results are presented in Table |[ 
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Figure 7. Shifts in the mean and fitted W mass as functions of the inverse weight 
cut off, 1/Vmax, for weight schemes Vi, V3 and V5, and y/s = 175 GeV. 



As mentioned above, for computational reasons, events with very large weights 
have been discarded. The dependence on the cutoff value, Kiiaxi was eliminated by 
calculating the multiplicity and mass shifts for three values of l^nax (20, 40 and 80) and 
then extrapolating to infinite cutoff. This method seems to be more reliable and less 
vulnerable to fluctuations than direct calculation with very high Knax- Fig. Q shows 
the values of the mean W mass, M^''^'', and the fitted M^* as functions of 1/T4iax. 
For Vi, the extrapolated value depended on the specific way the extrapolation was 
performed, and this ambiguity was added to the error shown in Table 4. 

From these numbers one can draw the following conclusions: 

• There is a clear correlation between the BE-induced shifts in the W mass and in 
the charged particle multiplicity in W+W" production: the larger the increase 
in charged multiplicity, the larger are the expected mass shifts. 

• Both V3 and V5 weights result in fairly small mass shifts. They are well-behaved 
numerically and probably give quite reliable estimates of the effect. The spread 
of values using various weighting schemes can be considered as indicative of the 
systematic errors inherent to this approach. 

• The fitted value for the W mass is less sensitive to BE effects than the mean 
over the full distribution, which has been used to estimate the effect in previous 
investigations Q. The estimated values for the shift in the fitted mass are less 
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Table 4. Values of differences in multiplicity and mass of the W boson for events 
with and without interconnecting Bose-Einstein correlations between the two Ws. 





Vi 


V3 




175 GeV 


3.8 ± 0.5 


1.8 ± 0.2 


1.0 ±0.2 


192 GeV 


3.7 ± 0.5 


1.7± 0.2 


0.6 ±0.2 


^j^^meon (MeV) 

175 GeV 


75 ± 15 


22 ± 11 


20± 14 


192 GeV 


92 ± 16 


34 ± 11 


38 ± 14 


AM^' (MeV) 
175 GeV 


12 ±9 


11 ± 7 


4± 12 


192 GeV 


15 ±8 


13 ± 7 


6±9 



than 20 MeV, implying that BE correlations are not too dangerous for the W 
mass measurements at the expected level of accuracy at lep 2. For the shifts in 
the mean W mass, values of order of a few tens of MeV were found, of the same 
general magnitude as in |l4|. In all cases the shift is towards larger masses, as 
expected on general grounds [T^ . 
• For all weighting schemes, the shift in Afyv increases with energy in the energy 
range considered, but the increase is fairly small. 

It is interesting to compare these results to the predictions based on the 
implementation of Bose-Einstein effects by shifting the momenta of final state particles 
p^ . The most important difference is in the particle multiplicity: event weighting 
naturally leads to an increase of the average number of particles due to Bose-Einstein 
correlations, while the momentum-shifting method assumes that the multiplicity is 
unchanged. The energy dependence of W mass shift is also different. The strong 
energy dependence in momentum-shifting scheme is a combination of two effects: the 
increase of the systematic shift for low momentum particles in the direction of smaller 
W momenta, and the differences in momentum spectra of W decay products for various 
energies, as stressed in Q. This seems to be less pronounced in the present approach. 

The study |l| confirms that the systematic effect of BE correlations on the W 
mass determination can potentially be quite large, as found in [|lj], although the 
actual values of the mass shift found here are somewhat smaller. The size of the shift 
is however quite sensitive to the procedure used to extract the value of the W mass. 
In particular it is observed that a fit to the lineshape of the W mass distribution has a 
much smaller systematic error from Bose-Einstein correlations than the average mass, 
due to the fact that the main effect on Mw in our scheme arises from the tails of 
the mass distribution, which contain very little information about the peak position. 
Hence it seems possible to keep the systematic error from this source below about 20 
MeV. Careful work linked to the actual event selection and fitting procedures used by 
the LEP experiments is obviously needed in order to assess this in the individual cases 
and to optimize the analysis procedures. Since the value of the mass shift is always 
positive (as also expected on general grounds), a further reduction of the systematic 
error by a factor two is in principle possible by assigning the expected shift as a 
correction to Mw- 

The comparison between hadronic decays in the W^W^^qqqq and 
W+W^— >qq£F^ channels gives a unique possibility to investigate the influence of 
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Bose-Einstein correlations on various properties of final state particles, such as mul- 
tiplicity, transverse and longitudinal momentum spectra, resonance properties and 
reconstructed jet characteristics. It is possible, that by taking proper care in the fit- 
ting procedures used, one may at the same time be able to use a large part of the 
hadronic W"'"W~ events for the W mass determination, and to study the interconnec- 
tion effects in the relatively clean setting of e+e~ W+W~ events, by restricting the 
study to the region of large, off-peak W masses where these effects are expected to be 
the largest. 
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2.3. Four- jet events^ 

The QCD processes e"'"e~ — s- (TP /j)* — > qqqq, qqgg form significant backgrounds 
( ^ 20%) to W'^W" — > 4 jet production at lep 2. It is important, both for the 
threshold and direct reconstruction Afyv measurements, that these backgrounds are 
well under control. In particular, the qcd Monte Carlos (mcs) used in the Afw analyses 
should correctly describe the relevant features (for example, the overall rate and the 
kinematic distributions) of the four-jet final states which pass the W"'"W~ selection 
criteria. 

In this connection it is worrying that certain aspects of four-jet production are not 
well described by the standard 'parton-shower (ps) -l-O(as)' mcs (jetset, herwig, 
. . . ). As discussed by G Cowan at this Workshop in particular Fig. 9, see also 
p3[), four-jet studies performed by the aleph collaboration at lep i reveal significant 
disagreement between data and mcs for distributions in the standard four-jet shape 



variables xbz, 4>ksw, ^'nr ^^'^ "^34 (f^'' definitions, see for example Ref. |23|). This 
suggests that the mcs do not provide a correct description of the angular correlations 
between the quark and gluon jets. On the other hand, 0{al) matrix element models 
(e.g. the JETSET 0{al) + string fragmentation model) give a much better description 
of four-jet final states The problem here is that matrix element models with 
'added-on' hadronisation cannot be reliably extrapolated from lep i to lep 2 energies 
- the hadronisation tuning is only valid at the lower energy. (It was reported at the 
Workshop that aleph have a special '©(a^) + PS + cluster hadronisation' version of 
HERWIG , but this is so far not publicly available.) 

Two studies directly addressed this problem at the Workshop. The first 
investigated whether qcd events which pass the W+W~ event selection do in fact 
populate the ranges of the four-jet angular variables where the mcs are known to have 
problems describing the lep i data. Fig. || shows simulations of QCD and W+W~ 
events at 172 GeV binned according to the four angular variables listed above. Two 
event selections have been used, giving similar results. The first, used for the workshop 
itself, was a linear discriminant constructed by aleph and designed for selection of 
totally hadronic decays of W+W" pairs at 161 and 172 GeV. Subsequently, this has 
been replaced by a more generic cuts-based selection; by energy-scaling the appropriate 
cuts, this has also been run on jetset ps Monte Carlo and real data at the Z° peak, 
and the discrepancies observed by aleph are seen to persist. The problem regions can 
be roughly characterised as follows: 

I cosxBzl, I cos^nrJ < 0.7, I cos (^Ksw I, I cos 0:34 1 < 0.5. (9) 

The plots shown in Fig. || correspond to generic cuts: the qcd predictions are taken 
from PYTHiA and the W+W~ predictions are obtained from koralw at 172 GeV. It 
is clear that both the selected W+W^ events and the predicted qcd background do 
populate the regions of concern. For example, approximately half the qcd events have 
I cos 034! < 0.5, a region where the ps mcs overestimate the lep 1 data by up to 15% 

The second study attempted to address the question of the origin of the 
disagreement between the ps mc and matrix element predictions for the angular 
variable distributions. Perturbative QCD predicts very specific angular correlations 
between the four final-state partons in c+e^ qqqq, qqgg (see for example Ref. p^). 

§ Prepared by R Jones, S Moretti and W J Stirling 
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Figure 8. Simulations of QCD and W+W events at 172 GeV binned according to 
the four angular variables discussed in the text 



These correlations are naturally included in a full matrix element calculation, but are 
not necessarily included in a PS emulation of the four-jet final state. 

To illustrate these correlations, we show in Fig. ^ the distributions in xbz (at lep i, 
using the Durham algorithm to define the four-jet sample) calculated for the qcd qqqq 
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Figure 9. Differential distributions in the Bongtsson-Zerwas angle 



(solid histogram, exact matrix element) and qqgg (dotted histogram, triple-gluon 
vertex graphs) processes. The former peaks around 90°, indicating the preference 
for the plane of the two secondary quark jets to be orthogonal to the plane of the 
two primary quarks. In contrast, the two secondary gluons prefer to be produced 
in the plane of the primary qq pair. Now at lep i the qqgg contribution dominates 
the total four-jet rate and so the | cosxbzI distribution is sharply peaked at 1. The 
maximum deviation between the data and the jetset mc occurs for cosxbz = 0, 
which is precisely the region where the qqqq contribution is maximal. It is possible, 
therefore, that the mc does not correctly include the qqqq angular correlations. To 
study this further, we can construct a toy matrix element ps-like calculation in which 
the correlations are switched off for the qqqq final state, i.e. the secondary g qq 
splitting is azimuthally symmetric about the gluon direction. The corresponding 
Xbz distribution is shown as the dashed histogram in Fig. ^. The distribution is 
significantly flatter, as expected. Also shown in the figure (dot-dashed histogram) 
is the prediction of a decorrelated version of the qqgg matrix element (again only 
including the triple-gluon- vertex diagrams). Since evidently the mc gives a good 
description of the lep i data for | cos xbz | '--^ 1 |^ we may conclude that the qqgg 
correlations are correctly implemented. (This is not perhaps surprising, since the 
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Figure 10. The ratio in Eq. ^) 



'abelianized' qqgg predictions do not differ markedly from the qcd predictions |2£ 
Fig. nn shows the ratio 



R = 



rfcr(qqgg) + dcr (exact qqqq) 
(icr(qqgg) + (iCT(decorrelated qqqq) 



(10) 



as a function of | cosxbzI (soHd hne). The prediction has the same qualitative features 
as the data/MC ratio of the aleph analysis This suggests that the lack of correct 
angular correlations in the qqqq part of the PS mcs is at least partly responsible for 
the disagreement with the lep i four-jet data. However the difference between data 
and MC seen by aleph is quantitatively larger than the ratio shown in Fig. |l^ ||^, |2^ . 
In fact it appears to be similar in magnitude to what would be obtained by either 
switching off the qqqq contribution entirely (dashed line in Fig. or giving it the 
same xbz dependence as the qqgg contribution (dotted hne). Further details of this 
study will be presented elsewhere [EgI. 
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3. Experimental issues in the measurement of M^y 

In this section various aspects of the experimental techniques and problems in the 
determination of Mw are discussed. For the 1996 lep data two methods were 
used. These are the measurement of the W+W~ cross-section near threshold {y/s = 
161.3 GeV) and the direct reconstruction technique for data at y/s = 172 GeV. These 
two methods are considered in turn. The list of topics addressed at the workshop is 
by no means comprehensive. 



3.1. Threshold method^ 



The main issue addressed during the workshop was the influence of interference 
between W+W^ production and other four-fermion final states. This is discussed 
in sections 3.1.1-3.1.3. In section [3.1.4 the dominant systematic uncertainty resulting 
from QCD background, is reviewed. 



3.1.1. Introduction The measurement of the W+W" cross-section at lep is the 
measurement of the cross-section for the process : e+e" W+W" fif2f3f4. This 
involves the production of two resonant W bosons (via the so-called "CC03 diagrams" 
- see Fig. 11). However, identical final states can be produced through different 
intermediate states: singly resonant W production (Fig. ^(a)), neutral current 
diagrams (Fig. |l2|(b)) and diagrams containing ^-channel W boson exchange (Fig. p^ ), 
can all contribute (see Note that the set of diagrams contributing to the muon 
semi-leptonic final state and the tau semi-leptonic channel are identical. These extra 
diagrams contribute not only as a background but, as the final state is identical to 
that obtained through the CC03 diagrams, there is also interference between these 
processes: one must sum the matrix element amplitudes, not just consider the squares 
of the amplitudes. 

The WW analyses are aimed at selecting events of the type resulting from the 
CC03 diagrams; hence the lep experiments generally choose to interpret their results 
in terms of the CC03 cross-section, rather than the full four-fermion cross-section. 
A further significant disadvantage in interpreting the results in terms of the full 
four-fermion cross-section is that this quantity is divergent in several of the Monte- 
Carlo generators used by the analyses (due to their massless fermion treatment). 
Although the CC03 cross-section is not strictly a gauge invariant quantity (see p7t), 
in practice it is well-defined and allows comparison between the measurements made 
by the experiments. Clearly, only the doubly resonant CC03 subset of diagrams have a 
production cross-section sensitive to the W mass around the threshold energy. Hence, 
one does not lose sensitivity in obtaining the W mass from the CC03 cross-section, 
rather than from a four-fermion treatment. 

In this section the methods adopted by the four lep experiments to compensate for 
the interference effects are considered, and a quantitative comparison of their results 
is attempted. 



3.1.2. Correction method All the experiments used four-fermion generators to 
produce simulation samples for the set of all diagrams and the CC03 subset. These 
cross-sections are defined as aaii and (Jqcq^ respectively and are used to make a 

§ Prepared by C J Parkes, P B Renton and M F Watson 
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Figure 11. Doubly resonant CC03 set of diagrams. 




(a) (b) 

Figure 12. (a) Singly resonant W production, (b) Additional diagram for producing 
particle anti-particle pairs 

correction from the measured (four-fermion) cross-section to the CC03 cross-section. 
In addition Delphi have used samples for the non-CC03 four-fermion diagrams {a4fbck)- 
The selected cross-sections after experimental cuts are denoted by cr' = ecr, where e is 
the efficiency. The following methods were used by the four lep experiments in their 
161 GeV cross-section papers |9[ where details of the procedures used can 

be found. The aim of all these procedures is to obtain the CC03 cross-section. Here we 
give the form of the correction term applied by each experiment: i.e. the term which 
contains the interference effect, as assessed from Monte-Carlo simulation. 

ALEPH and OPAL applied additive corrections to the measured cross-section. The 
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Figure 13. Example t-channel gauge boson graphs for the ei/cQq' channel 



background from all four-fcrmion processes is defined to be a'^^ — CqcqS' 
be subtracted from the observed cross-section, f The background includes a 'virtual' 
contribution due to the interference between identical final states; in the eT7cqq' channel 
the interference is large and negative, so the ultimate sign of the correction could be 
either positive or negative: in fact in this case the effective background cross-section 
is negative. 

DELPHI and L3 applied multiplicative correction factors. For Delphi this is simply 
the ratio of the result that is obtained when neglecting the interference between 
diagrams to the desired CC03 final result, i.e. "'^(^Qg/cr^jj — ct^^jj^,^. For L3 the correction 
factor includes the total contribution of the four-fermion diagrams and is defined as 



3.1.3. Correction results The Delphi correction factors, as a function of -^/i, are 
given in Fig. It can be seen that in the /iF^qq' and qq'qq' decay channels no 
significant effects are observed; whereas there are significant effects in the eFeqq' and 
tUfivi channels near threshold. The results were found to remain stable within the 
statistical errors over a wide range of different event selections. Thus at 161 GeV 
DELPHI finds strong negative interference in the eFeqq' and tUiivi channels which 
corresponds to a shift of about 30 MeV in the W mass. 

A quantitative comparison of the four-fermion corrections applied by the four 
experiments is complex due to the differences in procedure. Using information from 
30| , p^ , and other details made available during the workshop, an attempt 
is made to put the results of the other experiments in a form similar to the published 
ALEPH numbers. For the results presented in Table ^ the opal and Delphi values have 
therefore been scaled by the inverse of the selection efficiency (ecCOs)"^- There was 
not sufficient information available to transform the L3 result to this form. The Delphi 
and OPAL numbers include only the four-fermion backgrounds to the same final state as 
is under study, whereas the aleph number includes in the correction all four-fermion 
final states. Whereas the other experiments choose to include the background from 
a Z and radiated gluons in their two fermion production, the aleph numbers include 
this in the four-fermion correction: this explains the large discrepancy seen in the fully 

f The corrections quoted by ALEPH have been divided by the the CC03 selection efficiency (ecCOs)"^ 
to produce a total, rather than a visible, correction cross-section. 



refs p8, 29, 
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Figure 14. DELPHI interference correction factors as a function of the centre of 
mass energy 



hadronic decay channel. Although the Delphi results show the interference effects to 
be relatively insensitive to the event selection, when one adds in the full four-fermion 
backgrounds this will no longer be true. Hence the results quoted in Table ^ will also 
be sensitive to the analysis cuts. 

In general, satisfactory agreement is obtained between the experiments. However, 
there could be a possible difference in the results of the two four-fermion generators 
used by opal . It is hoped that this study will stimulate further co-operation between 
the experiments on this matter. 



3.1.4- Systematics The four experiments agree on the source of the dominant 
systematic error on the WW cross-section: the precision of the simulation of the Z{j) 
hadronic decay events, qcd at lep 2 is considered elsewhere in this report ([^ and 
section 2.3; here the techniques used to assess this dominant systematic uncertainty 
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Table 5. Correction factors cr^jj — ctqqoS' rescaled by the CC03 selection efficiency. 
The full details of the comparison are given in the text. The fi7fqq' column is the 
effective average of the three leptonic channels. Note that the result marked f is 
not comparable with the other results in that column (as it includes an additional 
source of background). Numerical results of the correction factors for L3 were not 
available, but the qualitative pattern of the corrections is the same as for the other 
experiments. 



Experiment 


Generator 




Decay Channel Correction (fb) 








qq'qq' 




ei/eqq' 


iuiqq' 




ALEPH 


KORALW 


+140 t 






-51±15 


+14 


DELPHI 


EXCALIBUR 


+40±21 


-3±7 


-34±7 


-44±15 


+1+11 


OPAL 


grc4f 


+56±34 


+6±17 


-20±18 


-26±32 


-15+18 




EXCALIBUR 


+28±23 


+18±11 


-49±12 


-18±21 


-21+10 



are reviewed. 

The background systematics were assessed by the experiments from a comparison 
of events produced by different generators/fragmentation packages; for example opal 
compared events produced using the jetset, herwig and ariadne models. Simulated 
events were then compared with data, both at 161 GeV, and because of the limited 
statistics available, with data taken at other energies: for example Delphi rescaled 
its selection cuts and compared data and simulation at 130-136 GeV and at 91 GeV. 
The errors ascribed by the experiments are given in Table The values of the errors 
assigned are reasonably compatible. The L3 value of this error cannot be compared 
directly to the others, but it would appear to be compatible. 



Table 6. Systematic error on the threshold cross-section from QCD events 


Experiment 


Error 


ALEPH 

DELPHI 

OPAL 

L3 


5 % on background cross-section (data/simulation comparison) 
6.5 % on background cross-section (different generators) 
10 % on background cross-section 
11.3 % on background cross-section 
4 % on total cross-section 
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3.2. Direct reconstruction 

A number of topics were discussed at the workshop, and are reviewed below. In 
sections 3.2.1 -^^2^3 various aspects of kinematic fitting and mass estimation and their 
impact in the W mass resolution are discussed. In sections 3.2.4 - 3.2.5 we describe 
some of the possible biases on Mw which can occur. Then, in sections 3.2.6 - 3.2.8 the 
question of four-fermion interference effects is discussed. 



3.2.1. Kinematic fit methods^ The resolution on the jet-jet mass is too poor to 
make a precision measurement of My^. Fortunately, several methods are available 
to improve the resolution. The methods adopted are either scaling techniques, in 
which the energy of the hadronic decay products of the Ws are rescaled to the beam 
energy, or constrained kinematic fitting techniques, in which overall energy- momentum 
conservation is imposed. These techniques rely essentially on the fact that the jet 
angles are measured with good precision, but the jet energies much less so. Both 
methods result in a considerable improvement in the jet energy resolution. 

In the absence of additional jets from qcd processes, the events WW — > qq'qq' 
give rise to a four-jet topology. The effect of extra jets is discussed below. If no 
additional constraints are imposed, then the fit for WW — > jjjj (where j signifies a 
jet) is a 4-constraint (4-C) fit. The four fitted jets are combined into two jet-pairs, 
corresponding to the two W bosons in the event, and this results in two estimates of 
Mw per event. Three possible choices are possible for this pairing. The two measured 
masses from the 4-C fit are highly anti-correlated, and can then be combined into a 
single average mass using the error matrix from the fit. A second option is to use a 
5-C kinematic fit in which the two fitted masses for the chosen jet-pairing are required 
to be equal. This is somewhat unphysical, since the masses of the two W's will not, in 
general, be equal, but it is found that this procedure gives a good mass resolution. It 
is thus of interest to understand which of these choices gives the better final precision. 

In the analysis procedures the WW jjjj event is often first forced into a 4-jet 
topology. However, in many events there are additional jets due to gluon radiation. 
Treating all events as 4-jet events will result in some loss of precision on Mw- In the 
preliminary Delphi analysis events with 5 jets are also allowed. However, the pairing 
problem becomes even more severe, with ten possible combinations. 

The experiments have adopted different methods to tackle the pairing problem in 
the fully hadronic channel discussed above, aleph use a combination of the values 
from two different kinematic fits; a 5-C fit (Xsc) ^^'^ ^ which both masses were 
compared to a reference mass (Xref)- Using the combination xlc + '^Xref^ Monte 
Carlo studies show that the correct solution is chosen in 79% of the cases. 

In the DELPHI analysis all the solutions are retained and used to create an 
ideogram probability distribution as a function of Mw- Each solution has a Gaussian 
distribution, with a width equal to the corresponding 5-C fit error on the mass. 
The relative normalization of the solutions depends on the significance of the mass 
difference in the corresponding 4-C fit. Each event is assigned a purity P on the basis of 
the value of the quantity i^min^min, where £'min is the minimum jet energy and ^min the 
minimum jet opening angle. This quantity is chosen as it has a rather different shape 
for the WW signal events and the Z/7 4-jet background events. The background 
in the event is 1 — P, and contains no information on Mw- The most favoured of 
the three solutions is normalized to P. For each solution a likelihood distribution, 
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as a function of Afw, is obtained by convolving a relativistic Breit-Wigner with a 
Gaussian distribution for the solution in question, with a width obtained from the 
fitted error for that solution. The likelihood distribution for the event is the sum of 
these distributions, normalized as described above. For the events classified as 5-jets, 
all 10 combinations are considered in this way. The probability distributions for each 
event are then multiplied together to give an overall likelihood from which the best 
estimate of Mw and its error are obtained. 

In the OPAL analysis the 5-C fit with the best xic probability is used. This is 
estimated to contain 68% of the correct parings. In addition, in certain cases, the fit 
with the second best xlc used. This is estimated to contain 25% of the correct 

parings. The statistical problem of retaining more than one solution per event must 
then of course be handled. 

The analysis of the semi-leptonic channel is considerably simpler than that of the 
fully hadronic channel. The hadronic decay products are usually forced to form two 
jets. For the WW =e,/x) channel there is an undetected neutrino, so if no 

additional constraints are imposed then the fit is a 1-C fit. Imposing equal masses for 
the leptonic and hadronic decay products thus gives a 2-C fit. For the WW jjri^ 
there are one or more additional undetected neutrinos, so that the r lepton energy is 
poorly known, but the decay products can be used as an estimate of the r direction. 
A rescaling technique can be applied to the fit results, or directly to the measured jet 
parameters. However, the final resolution on Afw is worse in this channel than for 
£ =e or /i. 

The main difference between the experiments in their preliminary analyses of the 
WW ii^i^ii =e,/^) channel is whether the individual event errors from the 2-C fit 
are used. These are used by Delphi, but not so far by aleph and opal. A significant 
difference between the WW ]]iv and WW jjjj channels is that the missing 
neutrino leads to marked differences in the event topology, and resultant kinematic 
precision, depending on its momentum. 

There is also a variation in the extracted precision coming from the choice of 
function used to fit the mass distribution. A study using Delphi Monte Carlo data 
has been carried out for the /^i7^qq' and eFcqq' channels and some of the results 
are shown in Table 0. In each case about 90 samples of Monte Carlo events, each 
corresponding to about 10 pb~^, were analysed. From the distribution of the fitted 
errors the mean and rms are determined. It can be seen that the use of the event by 
event errors gives an improvement of ~ 15% in the mean error. The rms spread is 
also improved. 

Table 7. Results on the expected error from fits to Afw using various functional 
forms and methods. The units are GeV. 



method 


mean 


rms 


single Gaussian 


0.598 


0.197 


single Breit-Wigner 


0.595 


0.215 


BW*Gaussian, a free 


0.610 


0.212 


BW*Gaussian, event errors 


0.503 


0.088 
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3.2.2. ALEPH rescaling technique^ Preliminary measurements of the W-boson mass 
from direct reconstruction have been made for the 1997 Winter Conferences by all 
experiments, and are discussed elsewhere in these proceedings 0|. All experiments 
employed a kinematic fit to reconstruct the event invariant masses. For the preliminary 
mass measurements using the 172 GeV data, Delphi, L3 and opal use 5-C kinematic 
fits to reconstruct the fully hadronic final state. The aleph Collaboration employs 
the kinematic fit in a somewhat different way (referred to here as 4C + Rescaling), 
and this is outlined in this section. 

In the ALEPH analysis of hadronic events from the 172 GeV data, events are forced 
into a 4 jet topology, using the Durham p scheme, and a 4-C kinematic fit is applied 
to the energy and direction of each jet constraining the total observed energy to be 
2£'bGam (i-e. ignoring initial state radiation (isr)). aleph Monte Carlo studies have 
shown that, on average, a 4-C fit reproduces the true jet energies more accurately 
than a 5-C fit with equal W masses. Then, to each of the three possible combinations 
of jet-jet pairs, the invariant masses, toi2,TO34, are rescaled to obtain a more precise 
dijet mass. The rescaled mass for jet-jet pair (1-1-2) is defined to be: 

TOf2 = mi2 X £;bcam/^12 £^bcam X \J I ^ PI2/ (H) 

and similarly for the other pair. The symbols E and p refer to the energy and 
momentum of the pair respectively. Errors coming from the loss of particles, or bad 
measurement of their momenta, are largely cancelled in the ratio P12/-E12; i-e. the W 
velocities are being used to obtain a reduction of the errors. A simple application of 
2-body kinematics gives: 

™f2 = ™12 

demonstrating that the rescaled masses are not improved estimates of their respective 
measured masses, but are now inter-related. It can be shown that the correlation is a 
function of the two true W masses in the event, anti-correlated by their measurement 
errors which are likely to be sensitively related to any cuts applied to the two measured 
di-jet mass distributions. 

A sample of sixty events is selected in the data and the best combination of jet-jet 
pairs, found as discussed above, is retained. The two separate rescaled jet-jet mass 
distributions are formed and each is fitted to a simple relativistic Breit-Wigner (bw) , 
from which M^j^, ^-nd their respective widths are extracted. Mw is obtained 

by averaging the two extracted masses taking into account the correlation expected 
from the Monte Carlo {p = -h0.32±0.10). 

Since the data sample is small, the statistical error is taken from the rms spread 
of the BW fitted masses, M^j^, obtained from 75 Monte Carlo samples each 

corresponding to 10.6 pb~^, where the input mass was 80.5 GeV. This gives: 

rms(Af^i) rms(M^2) ™s(Mi|,J rms(Af5C) 
463 MeV 483 MeV 420 MeV 506 MeV 

where is the average of and M^^2j taking into account their correlation. 

Thus the 4C -I- Rescaling procedure produces a smaller rms by 17% than the 5C fit. 
The corresponding average fit errors in the 4C -I- Rescaling and 5C analyses from the 
Monte Carlo samples are 320 MeV and 350 MeV respectively. The final calibration 
procedure increases these rms-based errors slightly: the corrected error for the 4C -I- 
Rescaling case is 450 MeV, which is still smaller than the 5C uncorrected fit result. 



'12 



'34 



4i?bcam£^12 



(12) 
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3.2.3. Further investigation of the rescaling technique The aleph rescaling method 
has been appHed to opal Monte Carlo data. Three samples with different input W 
masses were considered, 78.33 GeV, 80.33 GeV and 82.33 GeV. Each was generated 
at a centre-of-mass energy of 171 GeV. Samples with different masses were considered 
in order to investigate any possible phase-space effects. 

Three different fits to selected fully hadronic final states were considered, a 4-C fit 
where the average mass was used, a 5-C fit, and the aleph rescaling method described 
above. In each case the jet-pairings which most closely matched the parton level W 
pairings were used. No background was included. One thousand samples of 100 
fully simulated events were considered. The rms of the values for Af-w obtained from 
the different methods were taken as estimates of the error on the fitted mass. The 
correlation coefficient between the two masses from the rescaling method was found to 
be p ~ 0.45. Calibration curves relating the measured to fitted mass were determined 
(separately for each fit) from the three different Monte Carlo samples, assuming a 
linear relation. The errors on the fitted masses were then scaled appropriately using 
the calibration curve to give the error on the measured W mass. Table ^ shows the 
errors on Afw for the different methods and for the different Monte Carlo samples. 

Table 8. Errors obtained using the OPAL analysis procedure for different kinematical 
fit procedures, for three values of Mw The errors on Afw are in MeV and correspond 
to a sample of 100 events. 



type of fit 


Mw = 78.33 GeV 


A% = 80.33 GeV 


Mw = 82.33 GeV 


4-C 


369 


396 


447 


5-C 


354 


370 


339 


rescaled 


362 


378 


339 



For the opal Monte Carlo samples, the rescaled method performs almost as well as 
the 5-C fit. No improvement of ~ 20% in the resolution is obtained. It is interesting to 
note that when phase space effects are less important (i.e for the Mw =78.33 sample) 
the performances of three different kinematic fits are very similar. 

In summary, the preliminary aleph measurement of M-^ from the 172 GeV data 
used 4-C kinematic fitted masses rescaled to the beam energy. The opal study does 
not confirm the aleph finding that there is an improvement of ~ 20% in this rescaling 
technique compared to the 5-C fit methodf . 

3.2.4. Systematic shifts on the W mass obtained by direct reconstruction^ The 
measurement of the W mass by kinematic reconstruction of the invariant mass of 
the W decay products requires several steps. Each of these may be conceptually 
well defined, but requires some operational definition, an estimator, for which there 
could be various options. The choice is motivated mainly by achieving high signal 
efficiency and the best resolution for the reconstructed mass. However these estimators 
may be biased, introducing mass shifts which have to be evaluated and corrected for 
using Monte Carlo events. Such corrections can introduce model dependences and 
systematic errors. Various sources of mass shifts are discussed below and the results 
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t A similar study by C Parkes, performed after the workshop using the DELPHI Monte Carlo data, 
reached similar conclusions. 
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are summarised in Table These studies are carried out using the aleph selection 
and analysis procedures. 

The first step in an analysis of Mw is the event selection. W+W~events can 
be selected using a multivariate analysis (e.g. a Neural Network) which gives the 
highest signal to background discrimination. An upper limit for AMw of 80 MeV was 
calculated, for a given selection procedure, as the difference of the reconstructed mass 
obtained using the 'standard' cut on the selection variable, from the reconstructed 
mass obtained applying no such cut. In addition a 'standard' reconstruction analysis 
was applied to Monte Carlo events with and without ISR. The difference in the results 
obtained for the two cases was « 350 MeV. 

Jet clustering can be performed using different particle association criteria and 
combination schemes, e.g the Jade or the Durham algorithm and massive (E) or 
massless (P) schemes respectively. All types of algorithms introduce mass shifts. 
Different association algorithms introduce similar shifts, but the different combination 
schemes produce substantially different results. The value of AMw quoted here is 
the difference between the results obtained using a P or an i? scheme. P type 
schemes associate a larger fraction of particles to the correct jet, hence also better 
reproducing the parent-quark direction and energy. E schemes, which assign masses 
to the particles, give, on average, shifts from the parent parton values in energies and 
angles which are larger than for the P schemes. 

The next step of the analysis is forming di-jet pairs. The pairing can be performed 
using information related to the di-jet mass itself; such methods give the highest 
efficiency for correct pair association but introduce biases to the measurement. The 
quoted value of AMw is the maximum difference between results obtained using 
pairing algorithms based on several different variables, e.g. angular separation of the 
jets, of a kinematic fit, jet charge, etc. 

In order to improve the mass resolution a kinematic fit can be applied, using the 
constraints of energy and momentum conservation (4-C fit). In addition, the exact (or 
Gaussian) equality of the two W masses of the event can be imposed (5-C fit). The 
difference in the reconstructed mass using the two types of fit is given in Table |^. 

Finally, to extract Mw from the reconstructed invariant mass, the signal is fitted 
with a function, such as a Breit-Wigner. However there is no unique form for such a 
function and, depending on the choice, different Mw and errors (from the fit) on Afw 
can be obtained. Functional forms used for the fit include convolving the Breit-Wigner 
with a Gaussian and including a phase space correction factor. The value AAfw in 
Table jo] is the maximum difference between results obtained using different functions 
and fit ranges; the maximum difference arises from performing the fit with or without 
the phase space factor. The shifts at generator level from using a phase space factor, 
or not, are about 150 MeV. At the detector level there is additional distortion of the 
line shape from the kinematic fit. 

As we can see from the Tabic ||, the mass shifts and hence potential biases can 
be large, but it is the uncertainty on the necessary corrections that is important for 
the final measurement. Hence, at any stage, it is preferable to choose an algorithm 
that is not only efficient but also as unbiased as possible, so that any associated errors 
are small. Often there are physics grounds on which the choice for the most unbiased 
method can be based. For example, for the jet algorithm one solution is to assign 
particles to jets using a P scheme and then recompute the mass of the jet using an 
E scheme. Since nearly all the corrections rely on how well the Monte Carlo events 
describe the data, it is clearly important that this agreement is tested as extensively 



33 



Table 9. Breakdown of possible bias components. AMw is the difference between 
the results for Mw obtained with different analyses (not the difference of the result 
obtained from the generated M^) A 'standard' analysis was performed changing 
every time only the relevant component. 



Component 


Evaluation Method 


AA/w 


Event selection 


Change cut value 


< 80 MeV 


ISR 


off-^on 


350 MeV 


Jet algorithm 


Massive— >Massless scheme 


350 MeV 


Jet pairing 


Change algorithm 


Ri 100 MeV 


Kinematic fit 


4-C^5-C 


Si 300 MeV 


BW fit 


Change functional form 


< 300 MeV 



as possible. 

3.2.5. Linearity and biases The methods of extracting Mw so far adopted rely 
heavily on calibration from Monte Carlo generated samples which are analysed with 
the same procedures as the experimental data. In order to ensure that this calibration 
does not introduce a bias, Monte Carlo samples for a range of Mw values are used. 
The relationship between the fitted mass value Mf^^ and Mtrue is determined. This is 
found to be essentially linear in the region of interest and can be cast in the form 

Mfit - 80.35 = a(Mtruc - 80.35) + h . (13) 

Ideally one would like the linearity a — \ and the bias h — Q. Table shows the values 
of a and h obtained for the preliminary lep results at the 1997 Winter Conferences. 

Table 10. Available values of the linearity a and bias b for the different channels. 



channel 


cxpt. 


a 


b 


qq'qq' 


ALEPH 


0.938±0.035 


0.19±0.02 




DELPHI 


0.98±0.07 


-0.10±0.04 




OPAL 


0.953±0.012 


0.095±0.014 




OPAL(no ISR) 


0.975±0.007 


-0.077±0.008 


/iF^qq' 


ALEPH 


0.894±0.048 


-0.25±0.03 




DELPHI 


0.80±0.06 


-0.25±0.03 




OPAL 


0.89±0.01 


0.37±0.01 


eFeqq' 


DELPHI 


0.84±0.10 


-0.43±0.05 




OPAL 


0.905±0.009 


0.24±0.01 


£17^ qq' 


OPAL 


0.907±0.007 


0.33±0.01 




OPAL(no ISR) 


0.940±0.005 


0.07±0.01 



A value of a less than unity results in a corresponding loss of precision in converting 
Mfit to Mtruo- From Table |l^ it can be seen that this is most often the case. The opal 
values with and without iSR show that iSR is part of the cause for a being less than 
unity. Further work by the experiments is clearly needed to optimize the sensitivity. 
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3.2.6. Possible distortion by non-CC03 graphs The presence of non-CC03 graphs wiU 
change the shape of the fermion-fermion invariant mass distribution in the region 
around Mw compared to the case of CC03 graphs only. Ahhough this can be taken 
into account imphcitly in the analysis by using a full four-fermion generator for the 
"calibration" procedure, it is of interest to understand the size of the effects and the 
agreement between different four-fermion generators. Firstly the shifts in a specific 
four-fermion generator, wphact, are considered. Then other studies, including more 
detailed detector and selection effects, are discussed. 

3.2.7. Generator level studies of distortion by noTi-CCOS graphs using wphact|| An 
investigation of the shifts between the value of Mw extracted from a full four-fermion 
treatment and the sub-set of CC03 diagrams has been carried out at generator level 
using wphactQ. In each case distributions were generated for Mw = 80.356 GeV for 
the semi-leptonic channels eFud and /uFud. The W width used was Fw = 2.098 GeV. 
The effects of collinear initial state radiation were included, but not final state 
radiation. Two types of mass distributions were analysed. First the "true" qq' 
invariant mass distribution, using the generated four momenta of these particles, was 
used. For the second distribution the three-momenta of the u,d and I were used to 
compute the missing momentum Pmiss- This was used as an estimator of the neutrino 
three-momentum. The neutrino energy is taken as the modulus of Pmiss- The invariant 
mass of the tU^ system was then calculated, and will be referred to as m{£i'),-cc. 

An example of a generated mass distribution for the el^ud final state is shown in 
Fig.|l^ for the full four-fermion case. Some loose event selection criteria have been 
imposed. The minimum of the u, d and £ energies must be greater than 10 GeV and 
the ud invariant mass greater than 40 GeV. The angle between the lepton and the 
beam must be between 10 and 170 degrees, and the £u and £d angles greater than 5 
degrees. 

Fits have been made over the mass range shown in the plot to determine the W 
mass and width. A relativistic Breit-Wigner with a running width was used. The 
overall for the fits was poor, and the errors were multiplied by a factor of two 
to obtain those quoted below. This gives an acceptable x^/d.o.f. of about unity for 
the "true" distributions, but is still poor (x^/d.o.f. up to 2) for the reconstructed 
distributions. The results are the fits are given in Table |l|. 

Table 11. Results of fits to -A/w, using the WPHACT generator, showing the 
difference between the full four-fermion and CC03 diagrams for the "true" qq' and 
iUi masses, and for the reconstructed £u£ invariant mass distribution (see text). The 
units are MeV. 









distribution 


cuud 


fj,i/ud 


qq' true 


-3±1 


Oil 


ffi-f true 


-|-41±1 


Oil 


iu( recon. 


-|-40±1 


Oil 



It can be seen from Table |ll| that the differences in the fitted Mw values between 
the full four-fermion and CC03 cases for the "true" qq' distributions are small. For the 
channel eUud the shift is about 3 MeV, whereas for /uFud it is even smaller. The shift 
with respect to the generated mass is about 25 MeV, and on the width the change is 
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Figure 15. Distributions of the "true" qq' and "reconstructed" iU^ invariant masses 
from the WPHACT four-fermion generator for the eVud final state. The definitions 
of the true qq' and m(^i^)rec are given in the text. 



4 MeV or less. The "true" iUf fitted masses show no difference in the /iFud channel, 
but a 40 MeV shift is seen in the eFud channel. When the fits are performed on the 
reconstructed fUg distributions a similar pattern is observed: a shift is still found in 
the eVud channel, and no mass difference occurs in the ^Fud channel. The fitted 
widths also change significantly (increasing by 0.2 GeV), indicating that the shape of 
the distribution is distorted by the reconstruction procedure. This is also apparent 
from the poorer of these fits. 



3.2.8. Other studies of distortion by non-CC03 graphs^ This section reports on 
studies performed by the opal and Delphi Collaborations comparing CC03 and full 
four-fermion simulated events. 

In the opal study simulation events corresponding to an integrated luminosity 
of around 5000 pb~^ were generated. The full opal detector simulation and 
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reconstruction methods were then apphed. Events which pass the usual selection 
criteria are then analysed in two complementary ways. 

The first method applied by opal uses the true generated W mass Mtruc for each 
selected event. Mtruc is calculated by appropriately pairing up the generated four 
final state fermions as if they came from W"'"W~. Two invariant masses are then 
formed, and the average is found. A fit is made to the distribution of Mtruc to a 
function BW(M) x p{M), where BW(M) is a relativistic Breit-Wigner, and p{M) is 
a phase-space factor equal to the cm. momentum of a pair of W's of average mass 
AI. Then a comparison between the fitted mass parameter in the bw for CC03 and 
four-fermion cases is made, and the results for (four-fermion — CC03) (in MeV) are 
shown in Table |l2|for the different channels and for two four-fermion generators; grc4f 

and EXCALIBUR. 

Table 12. Differences between the values of Aftruc between four-fermion and CC03, 
for different channels and two different four-fermion generators. The units are MeV. 



generator 


qq'qq' 


qq'euc 


qq'fj.iy^{qq'TUr) 


grc4f 


-28±15 


-|-26±23 


-23±18 


EXCALIBUR 


-4±13 


-|-58±19 


+3±15 



Despite the large statistics in these samples, the results are still somewhat 
inconclusive. Both generators seem to indicate a positive mass shift for qq'eFo 
compared to the other channels, at the level of about 50 MeV, albeit with large 
errors. These results confirm the naive expectation that the shift is larger in the 
qq'el/e channel than the (\c{^iv^ channel, as additional diagrams contribute in the 
electron channel (see section 3.1.1) However, there also seems to be an overall shift of 
the grc4f numbers with respect to the excalibur numbers, by about 25 MeV. This, 
of course, is not expected. 

In the second opal analysis the difference Mrec - Mtruc, where Mroc is the 
reconstructed mass, is formed. From a fit to the peak region of this difference the 
mean is found. This is then a measure of the bias introduced by the reconstruction. 
Next the difference in this bias between four-fermion and CC03 is found. This study 
is done for events with only a small amount of ISR, because ISR distorts the resolution 
peak. The results are given in Table ^ 

Table 13. Mean of the differences between the values of Afrcc - Aftruc between four- 
fermion and CC03, for different channels and two different four-fermion generators. 
The units are MeV. 



generator 


qq'qq' 


qq'ej/c 


qq'Ati/^(qq'r;/T) 


grc4f 


0±12 


-26±20 


-30±20 


EXCALIBUR 


+23±10 


-5±16 


+28±16 



It is again difficult to draw very clear conclusions from this set of numbers. Naively 
one would have expected these all to be close to zero; and indeed they may, but only 
two lie within 1 s.d. of zero. Again, the grc4f numbers are systematically shifted 
compared to excalibur. Whether these differences stem from the original generators 
or from the modifications needed to implement them in the experimental environment 
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(e.g. Coulomb corrections, qcd corrections to Fw, an angle cut on the electron, 
treatment of masses etc.) is still an open question. 

The DELPHI study was made using the excalibur four-fermion generator to 
produce CC03 and four-fermion event samples. The event selection was made at the 
generator level, where cuts were applied that mimic those of the collaboration's full 
analysis. The samples were then analysed in a similar manner to the first opal method. 
The results of this study, and a further similar study in the Delphi Collaboration Q , 
are compatible with the opal excalibur results given in table The shifts in 
the qq'qq' and /xF^qq' channels are found to be small ( < 10 MeV). For the el^eqq' 
channel the shifts are small ( ^ 20 MeV) when fitting the qq' invariant mass, but 
about 75 MeV when fitting the tPg invariant mass. Thus it would appear that it is 
mainly in the tP^ invariant mass distribution that a distortion of the Breit-Wigner 
shape by non-CC03 graphs occurs. 
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3.3. Error from the determination oj the lep energy^ 

At LEP 1 the absolute energy scale is determined by the very precise method of resonant 
depolarisation, which gives an accuracy of better than 1 MeV in the average circulating 
beam energy. This technique becomes increasingly difficult at higher beam energies, 
so that the accuracy at high energy is largely determined by the highest value of the 
energy at which resonant depolarisation can be made. 

Several attempts were made in 1996 to achieve transverse polarisation at as high 
an energy as possible. A measurable level of polarisation was obtained at 50 GeV on 
several occasions, but so far attempts at higher beam energies have not succeeded. 

The 1996 beam energy was obtained with a precision of about 30 MeV |Q. 
By far the largest component of the error came from extrapolating the polarisation 
measurements at lower energies (45 and 50 GeV) to the operating beam energies of 80.5 
and 86 GeV. The extrapolation errors were 24 and 29 MeV respectively. The other 
significant components in the error came from the relative fill to fill normalization (10 
and 5 MeV respectively) and the uncertainties in the modelling of the RF corrections 
(5 and 5 MeV respectively). The lep beam energy error is a common systematic 
error when combining results from the four lep experiments. However, because of the 
limited luminosity for the 1996 data, the lep energy error remained a small component 
of the overall Afw error. 

In order that the error from the lep beam energy remains a small component 
in future high luminosity running it is desirable that it can be reduced to 15 MeV 
or better||. In order to achieve this the extrapolation error must be reduced. Two 
approaches are being pursued by the lep Energy Working Group. The first is to 
understand better the relative energy scales at 45 and 50 GeV, where measurements 
have been made. The relative error (cte) for these two energies in 1996 was about 
4 MeV; and this leads to the extrapolation errors quoted above. A significant part of 
this error comes from the rather limited number of such polarisation measurements. 
In only one fill was the energy measured at both these energies, and in addition 
there were two fills where measurements at 50 GeV were made. More depolarisation 
measurements at these energies should, provided there are no surprises, allow a 
significant reduction in cte- 

The extrapolation with energy in 1996 was taken to be linear. If cte can be 
substantially reduced, then possible non-linear terms in the extrapolation could 
then constitute a sizeable part of the extrapolation error. Thus obtaining resonant 
depolarisation at energies higher than 50 GeV is a high priority. During the 
1996/97 shut-down improvements were made to lep which are designed to produce 
fiatter orbits, and thus improve the polarisation build up at higher energies. 
Polarisation measurements at 55 GeV, or possibly 60 GeV, would help considerably in 
understanding the extrapolation. If polarisation at higher energies cannot be achieved 
then there is a possibility to get a third point somewhat below 45 GeV. However this 
possibility is limited by the range of the nmr devices which measure the dipole field 
and are needed in the extrapolation and modelling procedure. 

In summary, provided that there are no hidden surprises, and that sufficient time 
is devoted to making resonant depolarisation measurements, then a lep beam energy 
error of 15 MeV should be achievable. 



§ Prenared by P B Renton 

II In Ifil a target error of 12 MeV was used 
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3.4- Error extrapolation^ 

The report of the Mw Working Group in the lep 2 Workshop YeUow Book (YB) 
Q contains estimates for the Hkely attainable precision on Mw from both the 
threshold and direct reconstruction method. These estimates contain the most 
accurate information available in late 1995; they were based on Monte Carlo simulation 
work and theoretical studies. However, now that a significant amount of lep 2 data 
has been collected and Mw measurements have been made using both methods (for 
a summary see [Q), it is worth reviewing these estimates. In the following sections 
we compare the precision currently achieved by the LEP collaborations with that 
predicted in the Yellow Book study, and use this comparison to refine the predictions 
for future attainable precision on Mw as a function of luminosity. 

3.4- 1- Threshold method The method of extracting Af^v from the W+W~ cross- 
section at threshold is discussed in In 1996 the four LEP experiments combined 
obtained a total luminosity of approximately 40 pb~^ at a centre of mass energy 
of 161 GeV. In Tables p^jl6| we consider the three decay channels {tPtlvi, qq'tui, 
qq'qq') and compare the Yellow Book predictions (labelled YB) with the published 
experimental results of the four lep collaborations (labelled achieved); in the light 
of experience we then provide new estimates that would be obtainable for future 
running. It is important to note that great effort has not yet been made in reducing 
the systematic errors, since for the current integrated luminosity the statistical errors 
are completely dominant: indeed currently the statistical errors from the number of 
simulation events produced can be a significant component of the systematic error. 

Table |l^ shows the signal efficiency, background cross section, and systematic 
error for the tUiiui decay channel. The systematic errors obtained by the four 
experiments already largely match the YB estimates. In addition, aleph and opal 
report combined efficiencies for this channel in excess of the 60% YB prediction, 
however the backgrounds are somewhat larger than anticipated. To extrapolate 
to higher luminosities, we assume an efficiency of e = 75% and a corresponding 
background estimate of 0.02 pb. We assume a (common) systematic error of 4%, 
a value already achieved by one experiment. Using these values we may estimate the 
fractional errors on ctww for different luminosities. Table |l^ compares the estimated 
('Oxford') and YB errors for £ ~ 10, 50, 100 pb~^ per experiment. Note that for this 
channel the Oxford and YB estimates are almost identical, the increase in background 
being compensated for by a slightly higher efficiency. 

Table 14. Signal efficiency, background cross section, and systematic error for the 
£vilv(^ channel. 





O-sig 


0"bkg 


e 


(A(T/(7)sys 


YB 


0.23 pb 


0.01 pb 


60% 


2.7% 


Achieved 




0.02 - 0.06 pb 


40 - 75% 


4 - 30% 


Projection 




~ 0.02 pb 


~ 75% 


^ 4% 










(mostly common) 



The results of the ^^qq' decay channel study are given in Table |T^. The most 
significant difference between the results achieved and the YB predictions is in the 
substantial improvement in efficiency in the qq' rFr channel. The YB study made 

§ Prepared by D G Charlton, P Dornan, R Edgecock, P B Renton, W J Stirling, M F Watson 
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no attempt to optimize this channel, and as a resuh the efhciency was only 5%. In 
contrast, the lep experiments have already achieved 40 — 50% efficiency, albeit with a 
significantly larger background. For the projection to higher luminosities we assume an 
overall efficiency of e = 75% and a corresponding background of 0.08 pb. We assume 
an overall systematic error of 3% for this channel, of which 2% is taken to be common. 
Comparing the fractional errors on the qq'-fi^i cross section listed in Table |l^, we see 
a significant improvement in the predicted overall precision with respect to the YB. 

Table 15. Signal efficiency, background cross section, and systematic error for the 
qq'ii'l channel. 





O-sig 


0"bkg 


e 


(A(7/(T)sys 


YB 


0.76 pb 


0.03 pb 


47% 


2.7% 


Achieved 




0.08 - 0.19 pb 


60 - 76% 


3 - 10% 


Projection 




~ 0.08 pb 


~ 75% 


~ 3% 










( ~ 2% common) 



It is more difficult to make a similar assessment for the expected precision for 
the qq'qq' channel. After a loose preselection to eliminate the major component 
of the qq'(7) background, the four experiments use a combination of multivariate 
and sophisticated cut based analyses to separate the W+W" signal from the QCD 
background. In contrast, the YB predictions were based on a simple cuts selection 
, yielding a signal efficiency of approximately 55% and a purity of approximately 
70%. Given the nature of the analyses performed, instead of quoting backgrounds and 
efficiencies, we provide a statistical sensitivity both for the present results and for our 
predictions. In order to extrapolate to higher luminosities we take the current aleph 
statistical sensitivity of ±28% for 10 pb^^ as typical, and assume a common systematic 
error of 5%, slightly larger than the YB estimate, see Table |l6|. The corresponding 

Table 16. Signal efficiency, background cross section, and statistical/systematic 
error for the qq'qq' channel. 





CTsig 




e 


(A(7/(T)stat 


{Aa/a)sys 


YB 

Achieved 
Projection 


0.94 pb 


0.39 pb 


55% 


28% (ALEPH) 
~ 28% 


4% 

~ 5% 
(common) 



fractional errors on crww are listed in Table |18|. The improved signal selection efficiency 
leads to an overall slight reduction in the cross section error compared to the YB 
estimates. As anticipated in similar errors on the cross-section are obtained in the 
qq'iPg and qq'qq' channels. 

In order to translate a cross section error into an error on Mw several additional 
(common) systematic uncertainties need to be taken into account. These are 
summarised in Table As discussed in detail in Section 3.3, it is estimated 

that the lep beam energy error (currently ±30 MeV) can ultimately be reduced 
to ±15 MeV. A value of ±12 MeV was assumed in the YB study The 
theoretical error on the W+W" cross section also contributes to the error on Myj: 
AAfw ~ 17 MeV x (Aa-/a-)thy. In the YB study a ±2% error was assumed, arising 
in part from (unknown) higher-order loop diagrams involving Higgs boson exchange. 
However calculations performed since the lep 2 workshop |35j have shown that the 
corrections are intrinsically small and decrease rapidly with increasing Mjj. Using 
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the range of values for obtained Afjj from direct searches and global electroweak fits 
essentially eliminates the error from this source. The remaining theoretical error on 
uww is therefore expected to come from higher-order finite QED corrections, and so 
it is reasonable to assume an error of ±1% for future projections. 

Table 17. Additional common systematic errors for the threshold mass 
measurement. Note that AMyf AE^^^^^ and AAf^ 17 MeV X {Aa/a)^^^y. 





A-Bbcam 




YB 

Achieved 
Projection 


12 MeV 
30 MeV 
15 MeV 


2% 34 MeV) 
2% H 34 MeV) 
1% H 17 MeV) 



Combining these two systematic errors with the cross section errors discussed above 
gives the estimated overall errors on Mw at £ = 10, 50, 100 pb~^ (per experiment) 
listed in Table |l^. Also shown are the YB estimates. (It is interesting to note that 
the YB estimate of ±220 MeV for 10 pb^^ is exactly the current error from the 
combined lep experiments.) The key point to note is that a significant improvement 
on the projected YB errors can be expected. In particular, a relatively modest 
£ = 50 pb^^ total luminosity per experiment could now be expected to yield an 
error of AAfw — 86 MeV from the threshold cross section measurement. 

Table 18. Anticipated and new projected fractional errors on (tww a^nd 
corresponding overall combined error on Mw for three different total luminosities. 



luminosity/expt. 


channel 


YB 


Oxford 


10 pb~i 




34% 


32% 






19% 


15% 




qq'qq' 


20% 


15% 




combined 


13% 


10% 




AMw (total) 


220 MeV 


170 MeV 


50 pb-l 


&7££u£ 


15% 


15% 




qq'Wi 


9% 


7% 




qq'qq' 
combined 


10% 
5.9% 


8% 
4.9% 




AMw (total) 


108 MeV 


86 MeV 


100 pb-i 




11% 


11% 




qq'iue 


6.5% 


5.1% 




qq'qq' 


7.3% 


6.7% 




combined 


4.4% 


3.8% 




AAfw (total) 


84 MeV 


69 MeV 



3.4-2. Direct reconstruction During the second period of running in 1996 the lep 
experiments collected a similar total integrated luminosity (approximately AQpb~^) at 
an energy of 172 GeV to that obtained at threshold. However the future accuracy 
of the direct reconstruction method for measuring Mw is more difficult to assess: the 
analyses are currently at a relatively preliminary stage, and these data are expected to 
comprise a relatively small proportion of those provided by the full lep 2 programme. 
Much more work directed at reducing the systematic errors will be performed as 
more data are collected. Nevertheless, based on existing experience and the Yellow 
Book studies it is possible to make some educated guesses as to the likely ultimate 
precision. As a benchmark, we take the current aleph and Delphi (see [H) statistical 
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and systematic errors, summarised in Table based on approximately £ = 10 pb^ 
per experiment at 172 GeV. Note that the systematic errors currently range from 
60 MeV to 180 MeV depending on the channel and the experiment. The major 
challenge is to reduce these to the 0(20—25 MeV) level anticipated in see elsewhere 
in this report. As discussed above, a lep beam energy uncertainty of 15 MeV appears 
to be achievable. 

Table 19. Current statistical and systematic errors (in MeV) on A/w by direct 
reconstruction from the ALEPH (A) and DELPHI (D) experiments. The statistical 



errors are 


estimated 


errors. 




channel 


expt. 


Stat. 


sys. 




A 


510 


60 




D 


500 


94 


qqqq 


A 


450 


180 




D 


450 


75 



Table ^ summarises the expected uncertainties at the three luminosities C = 
100,300,500 pb^^ (per experiment). Also shown are the YB estimates. Note that 
the dependence of AAfw on the beam energy in the 180 — 200 GeV range is expected 
to be weak |^ . The estimates for the qq'qq' channel in Table ^ do not include any 
contributions to the errors from colour reconnection or Bose Einstein correlations. 

Table 20. Expectations for AMw (in MeV) combining all four LEP experiments 
for three different total luminosities. Note that the qq'qq' channel and combined 
channel estimates do not include any additional errors from colour reconnection or 
Bose Einstein correlations. 



channel 


luminosity/expt. 


Stat. 


sys. 


LEP 


total 




YB 


qq'£F< 


100 pb"^ 


71 


21 


15 


76 








300 pb-l 


41 


21 


15 


49 








500 pb-l 


32 


21 


15 


41 




44 


qq'qq' 


100 pb"^ 


71 


24 


15 


77 








300 pb-i 


41 


24 


15 


50 








500 pb-l 


32 


24 


15 


43 

(ecol.rec.a 


3BE) 


45 


combined 


500 pb-^ 








35 

(ecol.rec.g 


3BE) 


34 



Hence the Mw results obtained on the 1996 data sample (using both methods) are 
broadly in agreement with the YB predictions. The opportunity is available at lep 2 
to obtain precision measurements of A/w from these two independent methods. We 
conclude that there is every reason to believe that the direct reconstruction method 
YB estimate of AAf^v < C(50) MeV can indeed be achieved. 

4. Summary 

After one year of lep 2 data taking, measurements of the W mass have been made 
using two different techniques: from the threshold cross-section and from direct 
reconstruction. The results from these methods are compatible and the errors are 
dominated by their statistical components. A common theme of the various studies 
performed at the workshop was to ascertain whether the systematic errors were well 
understood. That is, are they understood at a level where one can extrapolate 
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with confidence to the full lep 2 statistics and ensure that systematic eff'ects will 

not dominate? 

The effects of colour recombination and Bose-Einstein statistics on the 
reconstructed W mass in the fully hadronic channel were explored in some detail. 
These studies supersede those made for the CERN lep 2 Workshop and in both cases 
suggest that the error could well be smaller than previous studies have indicated; 
although more work is needed for a definitive conclusion on these topics. The 
importance of relating the predictions of the various models to measurable quantities, 
such as multiplicities, was stressed. 

The background from the qcd processes e+e" — > {Z^/'y)* — > qqqq, qqgg to 
W+W~ 4 jets is large and needs to be well understood. One particular problem 
is the extent to which the precise lep 1 studies can be extrapolated to lep 2 data. A 
study at the workshop ascertained that the relevant W+W~ events do populate the 
ranges of the four-jet angular variables for which there arc discrepancies between data 
and Monte Carlo at lep 1 . An attempt was also made to understand the differences 
between the ps mc and matrix element predictions for angular variables. Part of these 
differences was found to be due to the lack of correct angular correlations in the qqqq 
part of the PS mcs. 

The correction factors in going from the measured W+W~ four-fermion cross- 
section to the CC03 cross-section, as determined by the lep experiments, were 
compared. This was difficult because of the different ways that these were implemented 
in practice. General agreement in these calculations was obtained. 

A detailed investigation was made of many of the systematic effects which need to 
be understood to make a precise measurement of Mw by direct reconstruction. These 
included studies of the problems of finding the correct pairing in fully hadronic decays 
and of the benefits and disadvantages of the various strategies of kinematic fitting and 
rescaling techniques. It was found that for these and other studies (on biases, linearity 
of response, distortion of the lineshape by non-CC03 graphs etc.) there was sufficient 
understanding at present but that further studies were needed to ensure the full lep 
statistics can be exploited. The lep energy error, which is a common systematic to all 
measurements of Mw can probably be determined with suflacient accuracy, provided 
the much needed studies show no surprises. 

Finally, a comparison of what can be expected if there was to be a further threshold 
run, and a review of the expected precision of the direct reconstruction method, were 
reported. From experience of running at 161 GeV, it appears that the precision on 
Mw from the threshold cross section method estimated in the Yellow Book study can 
be improved upon. Overall, provided there are on-going studies on the systematics, 
the desired goal of AMw < 0(50) MeV, seems achievable. 
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